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SIMUIATED MJCLEAR BEkTim OF LIQUID HYDROGEN IN A PROPELLANT TANK 


by Sidney C. Huntley and James W. Gauntner 
Levis Research Center 


SUMMARY 

An experimental study vas made to simulate nuclear heating of liquid hy- 
drogen in a propellant tank by using an electrical immersion heater and radiant 
heaters. Flow tests vere made in a 125-gallon tank pressurized to 2 atmospheres 
vith a flow rate of about 0.04 pound per second. Test results shoved that the 
inoperative immersion heater did not alter the exit temperature history. Oper- 
ation of the immersion heater tended to result in a more completely mixed liq- 
uid than existed vith nuclear heating. A comparison of bottom (radiant) heat- 
ing vith nuclear heating shoved a nearly identical generalized exit temperature 
history vith a similar generalized heat input rate distribution and equal heat- 
ing parameters^ although the initial total heat input rates varied by a factor 
greater than 6. High bottom (or liquid soirrce) heating vas observed to con- 
sistently res\ilt in a liquid disturbance near the surface d\n*ing the transient 
development of a temperature gradient. 


INTRODUCTION 

Interest in propellant heating has increased in recent years because of 
the desire to achieve minimum rocket vehicle veight. A knovledge of the ef- 
fects of propellant heating on such vehicle components as the propellant tank^ 
the pressurization system^ and the pump is desirable for optimization of the 
vehicle veight. Propellant heating in most rocket vehicles consists of a net 
inf line of heat to the propellant from thermal radiation and aerodynamic heating 
of the tank vails. Additional heating of the propellant in a nuclear vehicle 
results from nuclear radiation absorbed both in the tank vails and in the pro- 
pellant itself. The propellant heating problem arises from uncertainties re- 
garding the liquid behavior vhen subjected to one or more of the aforementioned 
heat loads. A knovledge of the liqioid behavior is required in making a compro- 
mise betveen tank press-ure level necessary to provide the net positive suction 
head requirements of the pump and the extent of thermal protection from the ex- 
ternal heating loads. 

In the early phases of the nuclear rocket program there vas little infor- 
mation describing the effect of vail and nuclear (or source) heating on liquid 
behavior. One of the first analytical studies of the problem, appearing in 
reference 1, consisted of a numerical approach to predict transient temperatiore 



profiles in the fluid either for closed tanks or for tanks with flow. There 
were no experimental nuclear data available for comparison* A qualitative in- 
vestigation of noncryogenic liquid behavior in a vented two-dimensional glass 
tank subjected to wall and source heating was conducted at MSA Lewis Research 
Center (ref. 2). An extension of this program consisting of an analytical and 
experimental investigation was made of the effects of nuclear heating on the 
temperature histories of liquid hydrogen. This investigation was conducted us- 
ing a three-dimensional tank representing a scaled-down version of a proposed 
nuclear rocket propellant tank (refs. 3 to 5) . Another experimental investiga- 
tion of side wall and bottom heating of liquid hydrogen contained in a tank 
similar to that of reference 3 was made and is presented in reference 6. The 
results of the work reported in references 2 to 6 show that the liquid behavior 
associated with side wall and nuclear (or bottom) heating tends to form a 
stratified region of warm fluid above a region of turbulently mixed fluid hav- 
ing a uniform temperature. An analysis was developed (ref ^ 5) and foxmd to 
predict satisfactorily temperature histories for cases with a major portion of 
the heating from below. The stratified region^ formed by the accumulation of 
boundary layer flow up the tank walls ^ was found to be dependent on the ratio 
of side wall heat flux to bottom heat flux. Low values of this ratio were 
found to yield thin stratified regions ; high values resulted in the entire liq- 
uid height being stratified. The effect of variations in this ratio on the 
temperature history of the liquid at the tank exit during constant press-ure 
outflow was also shown. In addition^ it was shown that the effects of total 
heating rate^ flow rate^ and tank pressure on the exit temperature history 
could be generalized. It was fiorther indicated that similar temperature his- 
tories were obtained from bottom heating as from nuclear heating provided that 
the heat input rate distributions were similar. 

The purpose of this report is to present the results obtained in an exper- 
imental investigation of electrical simulation of nuclear heating of liquid hy- 
drogen. The electrical simulation utilized an immersion heater to produce heat 
sources in the liquid from the heating of high-resistance 0.020-inch wire ele- 
ments arranged in a l/4-inch-wide by 1-inch-high grid network. Nuclear heating 
of the tank wall was simulated by using radiant heaters. This work was under- 
taken because of the complex and inflexible nature of performing tests in a nu- 
clear environment. The nuclear environment of reference 3^ the only available 
source of nuclear heating data with liquid hydrogen^ was selected as a specific 
case for simulation. Simulated nuclear heating tests were performed in the 
tank described in reference 6. Data were also obtained to determine the effect 
of the presence of the immersion heater (without power applied). The high bot- 
tom (radiant) heating work of reference 6 was extended to more nearly duplicate 
nuclear heating (without simulated source heating) and these results are also 
presented. The data are generalized in accordance with a set of scaling rules 
(presented in ref. 6) to make them comparable with other data obtained at dif- 
ferent operating conditions. 


APPARATUS 

The apparatus used in this investigation^ with the exception of the immer- 
sion heater^ was described in detail in reference 6; consequently^ only a brief 
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description is included herein. The test apparatus consisted of the test tank 
enclosed in a vacuum chamber with provisions for supplying liquid hydrogen^ 
pressurizing and venting the tank^ throttling the outflow of liquid^ providing 
source heating of the liquid with electrical elements^ and providing separately 
controlled heating of the tank wall and bottom (fig. l) . A vacuum system was 
provided which had the capability of producing a vacuum of less than 10“*^ torr^ 
thereby reducing the conductive heat leak through the vacuum space to a negli- 
gible amoimt. 


Tank with Radiant Heaters 

An exterior view of the tank with radiant heaters installed is shown in 
figure 2. The tank^ constructed of 304 stainless steely had a liq*uid capacity 
of about 125 gallons with a diameter of 32 inches. The tank geometry consisted 
of cylindrical walls with a 45^ half -angle conical bottom. Spherical sections 
were used to join the wall and bottom sections and to provide the bottom sec- 
tion with a small radius. A flange and transition section on top of the cyl- 
inder was provided for removal of the tank from the elliptical tank dome. The 
outer surface of the tank and the inner surfaces of both the wall and bottom 
heaters were sand blasted and spray painted with a flat black lacquer to pro- 
vide high emissivity s-urfaces for radiant heat exchange between the tank and 
heater surfaces . 

The radiant heaters were constructed from 304 stainless steel ribbons. 

The wall heater consisted of a double lead coil of 1- by l/l6-inch ribbon cy- 

lindrically wrapped with a 37— inch inside diameter and a 2— inch pitch. Each 
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coil had approximately 13 convolutions. The bottom ends of each coil were 
joined with a length of the heater material to form one continuous strip to 
which electrical power leads were attached at each end. The overall height of 
the wall heater was about 30 inches . 

The bottom heater consisted of a double lead coil of 3/4- by l/l6-inch 
ribbon helically wrapped with a 45^ half angle from the axis with a major diam- 
3 3 

eter of 36^ inches using a l“-inch pitch. Each coil had approximately 13 con- 
volutions with the bottom two lead ends joined to form one continuous strip. 

The double lead coils of one continuous strip were used for both the wall 
and bottom heaters to prevent inductive coupling between the heaters and the 
tank. Voltage -regiilated alternating current was separately supplied to each 
heater . 


Immersion Heater 

The immersion heater (fig. 3) consisted of 25 individually controlled 
electric horizontal heating elements^ each one being a continuous piece of 
0.020-inch wire^ having been wound on a form^ stress relieved^ and then epoxied 
to a frame of l/S-inch stainless steel wire which itself had been previously 
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stress relieved. A wire having essentially constant resistivity in the temper- 
ature range from 36^ to 500^ E was used. The 25 individual heater elements 
were assembled with the horizontal wires of each element at 90^ with respect to 
the horizontal wires of the two adjacent elements to deter the development of 
unnatural flow patterns in the liquid. 

A 1/2-inch clearance between the heater assembly and the tank wall was al- 
lowed to provide space for boundary layer movement and instrumentation leads. 
Because nuclear source heating yields a relatively smooth gradient of heating 
rate^ it can be more closely simiilated by heat sources whose horizontal and 
vertical spacings are very small; however^ practical considerations of con- 
struction led to the use of l/4-inch horizontal spacing between wires and 1-inch 
vertical spacing between elements. The nuclear source heating load (heating 
rate per unit volume) of reference 3 decreased exponentially with height of 
liquid^ and the heating load half way up the tank was only about 10 percent of 
the heating load at the bottom. Consequently^ the heater was only extended 

about half way (24-^ in.) up from the bottom of the tank. 

A constant heating load in the radial direction was assumed to be an ade- 
quate simulation of the nuclear heating load experienced in reference 3. 


Instrimientation 

Two temperature meas-urement systems were used for the test: one consisted 

of carbon resistors to measure temperatures in the cryogenic temperature range 
and the other of copper constantan thermocouples to measure gas temperatures. 
Thirty carbon-resistor thermometers were mounted on a rake extending down the 
tank axis. These thermometers^ spaced at l/4-inch intervals near the liquid 
surface of the full tank and at 4-inch intervals elsewhere^ were used to observe 
cryogenic temperatures. A thermometer placed at the tank exit was used to ob- 
serve the liquid temperature leaving the tank. Temperature data are presented 
in this report for thermometers located at 0^ 16.9^ 24.7^ 32.7^ and 40.7 inches 
above the exit port along the tank center axis. Additional data for one test 
are shown for thermometers spaced at l/4-inch intervals near the initial liquid 
surface. Thirty-six copper-constantan thermocouples^ mounted along the entire 
axis of the tank at 4-inch_^ or less^ intervals on the axial rake^ were used to 
observe gas temperatures. Three thermocouples were attached to the inner sur- 
face of each radiant heater to measure heater temperature. Previous experience 
with the carbon resistors gave indication that an accuracy within 0.1^ R was 
obtainable in the range of liquid hydrogen temperatures. The thermocouple ac- 
curacy was estimated to be within 1^^ E at temperature levels greater than 
138^ R. 

Other types of measurements included tank pressure^ liquid outflow rate^ 
and liquid level position. 


4 



PROCEDURE 


The desired temperatizres on the radiant heaters were established to obtain 
both the desired wall heat flux and the desired bottom heat flinc into the tank. 
These temperatures were adjusted by varying the electrical input power to each 
heater. If the particular test required soiirce heating^ power was applied to 
the immersion heater so as to achieve a specified gradient of source heating 
rate per unit volume as a function of distance from the bottom of the tank. 

This gradient was similar to the nuclear heating deposition gradient which was 
experienced in the nuclear tests (ref. 3). After both sets of heaters were op- 
erating^ the tank was closed and 5 seconds later the tank was pressurized with 
hydrogen gas to 1 atmosphere above the initial press-ure level. The liquid hy- 
drogen was initially at saturation temperature corresponding to the initial 
pressure level. Next, the fill valve was closed (fig. l) . After the tank had 
been closed for 20 seconds and a stable pressure achieved, the shutoff valve 
(fig. l) was opened and the outflow of hydrogen began at a rate of about 0.04 
pound per second controlled by the throttle valve. Additional press-urizing gas 
was used as needed to maintain a constant tank pressure during the entire flow 
run. 


Both the double wall of the outer shell and the tub -like enclosure around 
the shutoff valve (see fig. l) were filled with liquid nitrogen to ensure a 
constant reproducible temperature which would yield a low ambient heat leak. 
Uniformity in the initial conditions was obtained by opening the throttle valve 
with the tank shutoff valve closed and flowing liquid hydrogen from the supply 
dewar to chill the venturi and to establish the throttle valve setting to 
achieve the desired flow rate for approximately 15 minutes prior to start of 
outflow from the tank. 

Test conditions for the several experiments are presented in table I. 

Each test run is identified by run number, the type of heating applied, and the 
average values of wall heat flux and bottom heat flux are shown. Also shown 
are the initial height of liquid, the initial liquid temperature, the tank 
press-ure with resulting saturation temperature rise, and the flow rate. The 
other items will be discussed later. 


DATA ANALYSIS 

The heat input rate from the radiant heaters was based on a calculation 
involving the heat flux between the heaters and the wetted tank wall. The 
total heat input rate was based on additional calculations involving the power 
to the immersion heater and the time -averaged enthalpy change of the liquid 
flowing from the tank exit. The latter calculations required that the radiant 
heater calculation be normalized to obtain the total energy input to the liquid 
over the time to outflow. Justification for this normalization was based on 
no-flow and boiloff tests presented in reference 6, and the method is extended 
here to consider the immersion heater. 
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Heat Input Rate from Radiant Heaters 


Each heater was considered to consist of three zones ^ each having a uni- 
form temperature measured by a thermocouple in that zone. 

The net radiant heat exchange was then calculated between each zone of a 
heater and the corresponding tank surface^ assuming that both the heater and 
the wall were gray s-urfaces and the close spacing of the heater to the tank 
wall was analogous to concentric cylinders. The net heat exchange per unit 
tank area per unit time from each zone was computed from the eq^uation 

^ = ° f 

where a is the Stefan-Boltzmann constant and 



The symbols used in this report are defined in appendix A. Constant emissivi- 
ties of 0.8 were assimied for both surfaces. The heat input rate from the 
heaters was then obtained from 


Qh(^s) = 



T. cLS , 


( 1 ) 


where Xg is the height of liquid in the tank. Equation (l) will be normal- 
ized later as previously mentioned. 


Ihe time-averaged rate of heat input from the heaters was obtained by use 
of the expression 



(2) 


where tjn is the flow time reqmred to empty the tank and Xg(t) is obtained 
from the relation 


V(L) - Vt = 



( 3 ) 


where V is the average volumetric flow rate calculated from the slope of a 
straight-line curve fit of liquid volume history. Ihe liquid volume history 
was obtained by observing the time at which the liquid surface passed the sev 
eral temperature sensor locations. 
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Heat Input Rate from Immersion neater 


The total power delivered to all immersion heater elements at any given 


time (Q])j(t)) as measured by 




was assumed to be entirely trans- 


ferred to the lic[uid. As the liquid level passed below a given element, that 
element was disregarded in making the summation. The resistance of each ele- 
ment was measured prior to the run while the voltage supplied to each element 
was measured during each run. 


The average heat input rate from the immersion heater during the flow run 
was obtained from the expression 


*^m . 

Q]^(t) dt ( 4 ) 

The location and power history of each element along with the previous 
relation of liquid level with time (eq. (s)) gave the total immersion heater 

power as a function of liquid height (Q^^(xg)). 


Heating Rate of Liquid 

The average heating rate of the liquid during a flow test was based on the 
time -averaged enthalpy increase in liquid flowing from the tank exit- The av- 
erage increase in liquid temperature was first calculated from the expression 

1 / N 

at = / AT(0,t) dt 

tm Jq 

where AT(0,t) is the measured increase in exit temperature at time t after 
start of flow and tj^ is the time required to empty the tank. The average 
heating rate of the liquid was then obtained from the expression 

^ = Kp ^ ^ V (5) 

where p and c^ are the liquid density and specific heat at constant pres- 
sure evaluated at the average temperature + AT by using the National Bu- 

reau of Standards literatirre. The values of density and specific heat evalu- 
ated for each run are shown in table I. 

The heating rate distribution in the liquid was considered as being pro- 
vided by both immersion heater power and the heating rate from the radiant 
heaters. Because this assumption did not consider other heat inputs (such as 
radiant heat transfer from the dome, for instance), the radiant heater calcu- 
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lations were normalized to ottain equality among the time-averaged heating 

• • • 

rates (eqs. (2), (4), and (s)), that is to obtain Qz = + Qh- The heating 

rate of the liquid is consequently given by the expression 

• » 

• • - Q|vr • 

Qj(x3) = Qu(x3) + ^ (6) 

% 

Specific values of the average heating rate from radiant heaters Qg and the 

• • • 

heating error term (Q^ - Q[f)/Qjj for each run are shown in table I. 


Generalizing Parameters 


Generalizing parameters were derived in reference 6 to make test results 
applicable to a wide range of operating conditions. Ihe same parameters are 
derived herein, for completeness, along with an additional term to account for 
liquid heat sotirces. A system is assumed to consist of the liquid bounded by 
the tank walls and the liquid gas Interface at anytime during flow. The as- 
sumptions are made that there is no heat transfer across the Interface, that 
energy is time dependent and varies only in an axial direction, and that liquid 
properties are nonvariant and determined at a time -averaged temperature level. 
An energy balance of the system may be expressed in terms of the time rate of 
liquid enthalpy change, the rate at which enthalpy is being transported from 
the system and the total heating rate entering the system. In equation form, 
this is expressed as 



Q2(xg) 


( 7 ) 


vhere y is the -weight flov rate and V(xg) is the Yolume of liquid in the 
tank at any given time. 

A saturation heating rate Qg is nov defined as the energy required to 
heat the initial liquid content to saturation temperature in the same time re- 
quired to empty the tank. This satioration heating rate can he expressed as 


% = 


Kp cT V(L)AT, 


^m 


= K Cp w ATg 


■When equations (?) and (8) are combined, the energy equation then becomes 


*'v(Xg) ' ' 


Q. 




Q,(L) 


(8) 


( 9 ) 
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•where 


- 

ATg 

try, 

m 

H 

II 

II 

. 




( 10 ) 


In this form the energy equation states that the nondimensional tempera- 
ture history everywhere in the tank is a f-unction of the nondimensional heating 

rate distribution Q^(x )/Q^(L) and a heating parameter The heating rate 

distribution^ in some instances, originates from both heat transfer through the 
tank vail and from source heating in the liquid. "When both terms exist, it is 
convenient to express the initial magnitude of liquid source heating in terms 

of the initial total heating Qjj(l)/Q^(l) . Values of the initial total heating 
rate Q^(l)^ the heating parameter and the initial ratio of liquid source 
heating to total heating rate Q]^(L)/Q^(L) for each run are shovn in table I. 


RESULTS ARD DISCUSSION 

A series of tests vas performed to doc-ument the use of an electrical im- 
mersion heater as a means of simulating nuclear heating of liquid hydrogen. 

The experimental data obtained from this series of tests are presented. Test 
results are used to shov the effect of the presence of the inoperative immer- 
Sion heater on the liquid behavior. Test results of simulated nuclear source 
heating are compared with results from nuclear heating tests (ref. 3). Results 
of tvo tests without source heating but having a relatively high ratio of 
bottom-to-wall heat flux (hereafter called high bottom heating) are also com- 
pared with results from the nuclear heating tests . 


Experimental Data 

The heat input rate calculated from the heaters as a function of liquid 
height is presented in figure 4 for the several tests . The data points shov 
calculated values of heat input rate from the heaters obtained by using equa- 
tion (l). The heat input rate distribution for r-un 1 (fig. 4(a)) vas estab- 
lished by using a relatively lov bottom radiant heater temperat-ure and high 
vail heater temperature. This heat input configuration vas selected to compare 
with a typical high vail heating configuration with no immersion heater in- 
stalled (H - L configuration of ref. 6). The resultant average vail and bot- 
tom heat fluxes are given in table I. Relatively high bottom heating vas used 
for run 2 (fig. 4(b)); this configuration vas selected as representative of the 
type of heat input conf igioration which might be experienced in a nuclear rocket 
(see refs. 3 to 6) . R\m 3 (fig. 4(b)) is similar to run 2 in heat input con- 
figuration (within 10 percent) and is presented to illustrate the consistency 
of temperature history data. In the case of run 4, both the radiant heaters 
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and the immersion heater were used^ and the respective heat input rates are 
shown in figures 4(c) and (d). Also shown are the heat input rates for another 
test using both sets of heaters^ rim 5. For the latter two runs^ the magnitude 
of the heat input rates was about the same^ but for run 5 the immersion heater 
elements were selectively turned off as the liquid level approached each ele- 
ment. This was done to observe if a measurable difference in results existed 
when the heaters were left on or turned off and will be discussed subsequently. 

The temperature histories at several tank heights are presented in fig- 
ure 5. In general^ the temperature rise at any thermometer location increased 
after the start of flow. This increase in temperature rise with time shows an 
increase in temperature of the main bulk of liquid and^ as the liquid level ap- 
proaches a given location^ an increase in temperature due to a stratified layer 
near the surface. The high wall heat configuration^ rim 1^ shows a smooth in- 
crease in temperature rise at all locations (fig. 5(a)) with a trend toward an 
increasing temperature gradient with time throughout the height of liquid. For 
instance^ 400 seconds after the start of flow the temperature rise was about 
0.1^ F at the tank outlet but at a height of 16.9 inches the temperature rise 
was about 0.3^ F; at 24.7 inches^ 0.5*^ F; etc. For the high bottom heat con- 
figurations (runs 2 and 3; figs. 5(b) and (c)^ respectively) or the simulated 
nuclear heating configurations (runs 4 and 5; figs. 5(d) and (e)^ respectively)^ 
a near constant temperature region existed to near the liquid surface. A com- 
parison of figures 5(a) and (b)^ for instance^ shows that the region of turbu- 
lently mixed (constant temperature) liquid is more extensive with high bottom 
heating than with high wall heating. These results are typical of high wall or 
high bottom heat configurations and have been previously observed (ref. 6). 

It was seen^ with the present experimental data^ that the temperature at a 
thermometer location of 40.7 inches (initial surface level was about 44 in.) 
with high bottom heating had a tendency to increase. Then it levels for a 
short period (figs. 5(b) to (e)) shortly after the start of flow and be- 
fore the liquid surface dropped a significant amount. A similar tendency was 
not observed with high wall heating (run 1^ fig. 5(a)). A review of previous 
wall and bottom heating data (ref. 6) showed that this phenomenon was consist- 
ent. A discussion of the transient development of a temperature gradient near 
the surface is presented in appendix B along with additional temperature data 
near the surface during the first portion of test run 4. 


Effect of Immersion Heater Installation on Liquid Behavior 

Perhaps the most profound effect of the presence of the immersion heater 
on liquid behavior was expected at the liquid gas interface when the liquid 
surface passed each element. A disruption of the placid surface did occur dur- 
ing the flow of the liquid surface past the elements which could be seen from 
viewing the liquid with a television system. As the liquid surface drops past 
the element;, the surface is disrupted into small waves which disappear after a 
few seconds. This disturbance was observed to occur each time the surface 
passed a visible heater element. The immersion heater elements appear to have 
perturbed the liquid surface during outflow; however^ the significant effect of 
any change in liquid behavior should appear as a change in tank pressure or in 
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the exit temperatiire history. Evaporation and/or condensation vas anticipated 
due to the additional wetted surface of each heater element above that of the 
tank wall alone. However^ the liqiiid surface passing each element did not re- 
sult in a noticeable change in tank pressiore. 

A comparison of generalized test results with and without the immersion 
heater installed is presented in figure 6. Both the generalized heating rate 
distribution (fig. 6(a)) and the exit temperature history (fig. 6(b)) are shown 
for run 1 obtained with the immersion heater installed and for a similar run 
(H - L configuration from ref. 6) without the heater installed. A comparison 
of heating rate distributions indicates that the two cases were obtained with 
nearly the same generalized heat input. A comparison of the resultant general- 
ized exit temperature histories shows a general consistency in shape between 
the two curves^ which indicates the liquid behavior was not overly influenced 
by the presence of the immersion heater. The temperature history of run 1 

(^ = 0.444) shows a trend toward a higher temperature^ as expected^ than that 

• 

of the rim (ref. 6) without the heater (^ = 0.429). There was a noticeable in- 
crease in time during which the exit thermometer indicated saturation tempera- 
ture rise with the heater installed (fig. 6). This trend existed at other 
thermometer locations^ even at locations above the height of the heater (see 
fig. 5^ for instance). It is doubtful^ therefore^ that this difference re- 
sulted from the heater installation. 


Electrically Simulated Nuclear Heating 

The generalized test results with electrically simulated nuclear heating 
are presented in figure 7. The results obtained from test runs 4 and 5 are com- 
pared with the results of a nuclear heating test having a similar value of 
heating parameter (run 18.100^ ref. 3). Two heating rate distribution curves 
are shown in figure 7(a) for each run. The top curve depicts the total heating 
rate distribution both from heat transferred to the liquid from the tank walls 
and from heat sources in the liquid; the latter is shown separately by the bot- 
tom curve. The total heating rate of the liquid for run 4 (fig. 7(a)) closely 
matched that of the nuclear heating run although the amount of liquid source 
heating for run 4 was considerably less for two reasons: (l) the heating load 
(power per unit volume) applied to each element was about 40 percent less^ and 
( 2 ) the immersion heater extended up to about half of the initial liquid 
height- The results obtained from rim 4 (to be discussed later) made it unnec- 
essary to conduct further tests with increased heating loads. 

The low total heating rate of the liquid for run 5 (fig. 7(a)) at low liq- 
uid heights resulted from turning off the power to successive immersion heater 
elements before the liquid level came within several inches of any particular 
element. This was done in an effort to avoid perturbing the stratified layer 
near the liquid surface by power generation. The initial liquid source heating 
rate was about the same for both tests (runs 4 and 5)^ and consequently^ the 
difference in liquid height for any given liquid source heating rate indicates 
the approximate height of liquid above an element when its power was turned 
off. 
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The generalized exit temperature histories for the two simulated nuclear 
heating rims are shown in figure 7(h) and compared with the results of the nu- 
clear heating test. The results of all three tests appear, in general, to he 
about the same during the early portion of the tests. In the last quarter of 
the flow period, a more significant difference appears. The exit temperatures 
of the simulated nuclear tests show a tendency toward a more rapid increase in 
temperature. With nearly eqml total heating rate distrihutlosas between run 4 
and the nuclear data (see fig. 7(a)), consistent temperature histories were an- 
ticipated as a f -unction of heating parameter. That is, the higher the heating 
parameter, the greater the temperatvure rise at any given time from start of 

flow (see ref. 6, for instance). Although the heating parameter for the nude- 
• • 
ar data (^ = 0.457) was slightly higher than the present test run 4 (^ = 0.434), 
the early history indicates an exit temperat-ure rise which is less by a small 
amotint. This trend indicates more heat has been stored in the stratified layer 
during the nuclear rrm. A comparison of the exit temperature histories near 
the end of the tests shows the previoixsly stored heat in the nuclear case re- 
sulted in an increase in temperature rise greater than those of the present 
test. 


A comparison of the results of run 5 and the nuclear data shows similar 
trends as the previous comparison. In this instance, the difference in temper- 
ature histories occurred with about the same heating parameter but with differ- 
ent generalized heat input distributions. Apparently the heating rate at low 
liquid height for run 5, being less than the nuclear data, resulted in lower 
exit temperatvires near the end of the run. 

Q?he previous comparisons show that, in general, the nuclear heating could 
be simulated by the -use of an electrical immersion heater. The trend toward a 
more completely mixed liqiaid with the immersion heater was taken to indicate no 
f-urther purpose would be fulfilled by appl 3 d.ng a higher heating load to simu- 
late more closely the nuclear liquid source heating. An attempt to reduce the 
indicated mixing effect by turning off the power to each element before the 
siarface passed the element resulted in a change of heating rate distribution. 
Another feasible sim-ulation of the total nuclear heating rate distribution was 
indicated in reference 6 by the use of bottom heating. 


Comparison of Bottom and Nuclear Heating 

Two tests were performed (runs 2 and 3) in which the bottom and wall heat 
flux were established to approximate the heating rate distribution obtained 
with the nuclear heating data of reference 3. The resilLts of these tests are 
compared with results of a nuclear test in figure 8. A comparison of the gen- 
eralized heating rate distribution (fig. 8(a)) shows that the total heat rate 
distribution is about the same for all tests although in the nuclear case, 
about 70 percent of the heat input was deposited directly into the liquid. The 
initial total heating rate of the nuclear data was over six times greater than 
that of either of the present tests (see fig. 8(a)). A comparison of the gen- 
eralized exit temperatixre histories (fig. 8(b)) shows that the results of test 
run 2 and those of the nuclear test were almost identical. These three tests 
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also had similar heating parameters despite the factor of 6 variation in the 
total initial heating rates. Test run 3 yielded a generalized exit temperature 
history consistent with its 3-percent-lower heating parameter. 

A comparison of figures 7 and 8 indicates that bottom heating tests more 
closely simiilated the nuclear heating test resoolts than the electrical immer- 
sion heater tests. Apparently the liquid behavior associated with the bottom 
heating tests resulted in a turbulent flow field similar to that with nuclear 
heating. The high degree of attenuation associated with nuclear radiation in 
both liquid hydrogen and the tank material (along with the uniform nature of 
the nuclear heat deposition) appears to result in liquid behavior similar to 
that experienced with convective heat flow from the tank bottom. 

The major contribution to turbulence is perhaps more from bottom heating 
than the liquid heat sources in nuclear heating. The liquid behavior near the 
sixrface d-uring the early period of flow (see discussion in appendix B) sug- 
gested that high bottom heating consistently influences the establishment of 
the depth of the stratified layer. This condition also existed with about the 
same heat flvx on both wall and bottom (H - H configuration of ref. 6). It 
seems possible^ therefore^ that the nuclear heating of the tank material, as 
opposed to the liquid source heating, may have had an overriding influence on 
the establishment of flow behavior. 


SUMMARY OF RESULTS 

The following results were obtained from an experimental study to simxilate 
nuclear heating of liquid hydrogen in a tank by means of an electrical immer- 
sion heater and radiant heaters . 

1. The use of an electric immersion heater to simulate nuclear heat 
sources in the liquid tended to resiuLt in a more completely mixed state of liq- 
uid than with nuclear heating. 

2. Simulating the generalized heating rate distribution of nuclear heating 
by means of radiant heaters gave nearly identical generalized exit temperature 
histories with similar heating parameters despite lonsimilar operating condi- 
tions. 


3. The presence of the l/4- by 1-inch grid network of the immersion heater 
appeared to perturb the liquid surface while dropping past the heater elements 
but without affecting the exit temperatirre history. 

4. High bottom (or liquid source) heating was observed to consistently re- 
s-ult in a dist-urbance of liquid near the surface during the transient develop- 
ment of a temperature gradient. 


Lewis Research Center, 

National Aeronautics and Space Administration, 
Cleveland, Ohio, December 13, 1965. 
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APEEWDIX A 


SmBOLS 

Cp specific heat of liquid at constant pressure, Btu/(lh)(°R) 
E voltage to heater 
/ gray-hody shape factor 
K conversion factor, W/(Btu/sec) 

L initial liquid height in tank, ft 

Py tank pressure, psia 

Q heating rate, W 

SI ratio of heating rates, heating parameter 
q heat flux, Btu/(sec)(ft^) 

R heater element resistance, ohms 

r radius, ft 

O 

S tank surface area^ ft 

T temperature^ 

t time from start of flow^ sec 

V tank volume, ft 

V volumetric flow rate, ft^/sec 

V ratio of volumes 

w weight flow rate, Ih/sec 

X height from tank bottom, ft 

e emmissivity 

•S ratio of temperature rises 

p density, Ib/ft^ 

a Stefan-Boltzmann constant, Btu/(sec)(ft^)(°P^) 
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T ratio of tojnes 

Subscripts : 
b bottom 

H heater^ radiant 

i unit 

in initial 

Z liquid 

m maximum 

N liquid source 

s sat-uration^ sirrface 

T tank 

V wall 

Superscript : 


time -averaged value 



APPENDIX B 


TRANSIENT DEVELOPMENT OF TEMPERATURE GRADIENT 

In the discussion of figure 5 it was observed that the temperatvire indi- 
cated by the carbon-resistor thermometers near the initial surface level some- 
times had a tendency to increase and then level off for a brief period shortly 
after the start of flow. This tendency was observed only with heating rate 
distributions having high bottom heating rates. This phenomenon was explored 
more fully dxiring run 4 by recording each thermometer output at 0.054-second 
intervals from the time the tank was closed until several seconds after the 
start of flow. Twenty-five consecutive measurements of each thermometer were 
averaged to reduce the effect of 60 cycle noise in the signals resulting in av- 
eraged data points 1.35 seconds apart. 

The resultant temperature histories of several thermometers are presented 
in figure 9. Also shown is the saturation temperature history. Thermometers 1 
through 13 are spaced at l/4-inch axial increments; 14 is 4 inches below 13^ 
and 15 is 4 inches below 14 (ref. 6). Thermometer 1 is 43.69 inches above the 
tank bottom. The temperat-ures measured at thermometer 15 are typical of tem- 
peratures below this location. The tank was closed 20 seconds before the start 
of flow; pressurization started 5 seconds after the tank was closed. Closing 
the tank caused the tank to self pressurize until the addition of pressurant. 
Pressurization of the tank resulted in an increase in liquid temperature due, 
probably, to isentropic compression of liquid. A temperature gradient also de- 
veloped in the top 2 or 3 inches of liquid during the time interval from pres- 
surization to the start of flow. The initial liquid level (as indicated by 
saturation temperature after pressurization) appears to be slightly above the 
position of thermometer 1. After the start of flow, the general trend toward 
Increasing temperature with time at any given location is the resiilt of heat 
addition and a decrease in liquid level with outflow but disturbances seem to 
appear in the upper regions of liquid (above thermometer 14) causing variations 
in indicated temperature histories at certain locations. 

The first disturbance occurred about 2 seconds after the start of flow. 

The change in temperature level was experienced almost simultaneously throiagh- 
out the upper region and probably resulted from a slight adjustment of liquid 
volimie in the tank upon opening the shutoff valve and filling the outlet piping 
to the throttle valve. The volimie adjustment would be observable in the upper 
region by changes in temperature level because of the existing temperature gra- 
dient. No effect of volume adjustment is noticeable without the gradient. A 
second distvirbance occiurred about 12 seconds after the start of flow. This 
change in temperatirre level was again experienced simultaneously but only near 
the surface (thermometers 3 to 7). 

A third disturbance about 15 seconds after the start of flow appeared as 
temperature -level changes of longer diiration (until about 40 sec after start of 
flow) and progressed downward with time. The downward progression is observed 
by the initial temperature level change occurring at lower locations with in- 
creased time. The progression continues for a depth of several thermometer 
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locations in about the same interval of time required for the liquid surface to 
drop the l/i-inch distance between two thermometers. 

It is probable that the third distinrbance resulted from a quantity of waim 
liquid having traveled toward the surface (probably along the tank wall)^ t-urn- 
ing radially inward near the surface^ meeting with other warm liquid at the 
center^ and being forced downward before dissipating. This phenomenon was also 
observed in the two-dimensional tank tests of reference 2. Schlieren photo- 
graphs of the early history of flow tests (in ref. 2) showed a downward flow of 
warm liquid with either wall heating or a combination of wall and bottom (or 
nonuniform source) heating. However^ in the liquid hydrogen experiments (pres- 
ent data and ref. 6) the results of this downward progression were observed 
only with high bottom (or nonuniform source) heating. The transient develop- 
ment of a quasi-steady-state temperature gradient in liquid hydrogen is thus 
seen to be dependent on the heating rate distribution into the liquid. Wall 
heating of liquid hydrogen appears to contribute in a smooth manner to the de- 
velopment of temperature gradients which may extend over the major depth of 
liquid (see fig. 5(a)^ for instance). Bottom heating^ on the other hand^ con- 
tributes to a generally more turbulent flow of fluid which dist-orbs the wall 
heating contribution to the temperature gradient within a short time period 
(about 15 sec in the example shown) and resiilts in a thinner stratified layer 
above a region of turbulently mixed liquid. 
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II I 


TABLE I 

R\an 

Type of heating 

Average heat flux, Btu/(ft^) (sec) 
Wall, qL„ 

Bottom, 

Initial liquid height, L, in. 
Initial temperatures, T^^, °R 
Tank pressure, FT, psia 
Saturation temperatxure rise, 

ATg, ^ 

Average flow rate, w, Ib/sec 
Average density, p, Ib/ft^ 

Average specific heat, cT, 
Bbu/(lb)(OR) ^ 

Average heating rate from 

heaters, Qh, W 

• • • 

Heating error term, { 0,1 - Qw)/QH 
Initial total heating rate of 
liquid, Q^(L), W 

Ratio of.soirrce to total heating, 
^w(l)/Qz(l) 

Heating parameter^ ^ 


- TEST COroiilONS 


1 

2 

3 

4 

5 

Radiant 


Radiant and 




immersion 

0.0082 

0.0035 

0.0031 

0.0040 

0.0040 

0.0003 

0.0083 

0.0071 

0.0073 

0.0071 

43.94 

44.20 

43.49 

43.79 

42.31 

38.25 

38.25 

38.30 

38.20 

38.44 

33.78 

34.00 

34.22 

34.76 

34.67 

3.95 

4.00 

4.00 

4.22 

3.96 

0.0404 

0.0405 

0.0368 

0.0407 

0.0401 

4.310 

4.316 

4.297 

4.303 

4.296 

2.48 

2.54 

2.54 

2.57 

2.58 

63.15 

59.10 

52.34 

60.93 

64.94 

1.277 

1.654 

1.586 

1.486 

1.386 

184 

141 

121 

202 

200 

— 

— 

— 

0.308 

0.316 

0.444 

0.325 

0.306 

0.434 

0.462 
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Dimensionless temperature rise, i9 Dimensionless heating rate, Q^XjI/QjIL) 



(a) Heating rate distribution. 



(b) Exit temperature history. 

Figure 6. - Comparison of generalized test results with and without immersion heater in- 
stalled. 
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Dimensionless temperature rise, -3 Dimensionless heating rate. Q/cif(U 





Dimensionless temperature rise, -9 Dimensioniess 



(a) Heating rate distribution. 



(b) Exit temperature history. 

Figure 8. - Comparison of bottom and nuclear heating data. 
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Disturbance 
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Third 
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temper- 

ature 


Figure 9. - Temperature history during transient development of temperature gradient 
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aeronautical and space activities of the United States shall he 
conducted so as to contribute to the expansion of human knowl- 
edge of phenomena in the atmosphere and space. The Administration 
shall provide for the widest practicable and appropriate dissemination 
of information concerning its activities and the results thereof.^* 

— National Aeronautics and Space Act of 1958 
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